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ABSTRACT 

We derive [Cu/Fe] for 117 giant stars in ten globular clusters (M3, M4, M5, 
M10, M13, M15, M71, NGC 7006, NCG 288, and NGC 362) and find that glob- 
ular cluster Cu abundances appear to follow [Cu/Fe] trends found in the field. 
This result is interesting in light of recent work which indicates that the globular 
cluster iv Centauri shows no trend in [Cu/Fe] with [Fe/H] over the abundance 
range —2.0 < [Fe/H]< —0.8. Of particular interest are the two clusters M4 and 
M5. While at a similar metallicity ([Fe/H] ~ —1-2), they differ greatly in some 
elemental abundances: M4 is largely overabundant in Si, Ba, and La compared 
to M5. We find that it is also overabundant in Cu with respect to M5, though 
this overabundance is in accord with [Cu/Fe] ratios found in the field. 

Subject headings: Galaxy: abundances — Galaxy: halo — stars: abundances- 
stars: Population II — globular clusters: general 

1. Introduction 

Until a large-sample, high resolution study of field stars covering a wide range in 
metallicity revealed a definite trend in [Cu/Fe] with [Fe/H] (Sneden & Crocker 1988; 
Sneden et al. 1991a), the element copper had not received a great deal of attention in stellar 
abundance analyses. Sneden et al. found that [Cu/Fe] decreased sharply with [Fe/H], 
reaching [Cu/Fe] ~ — 1 at [Fe/H] ~ —3. Their reanalysis of other data in the literature 
along with averaged globular cluster abundances reinforced their field star findings. This 
work has been revisited by Mishenina et al. (2002), who confirmed this basic result with an 
expanded sample of halo field giants. 

Despite this clear trend in [Cu/Fe] with [Fe/H] among field stars, interpretation of its 
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origin is distinctly clouded. In the solar system, copper owes ~30% of its abundance to the 
s-process weak (massive star) and main (low mass star) components combined (Matteucci 
et al. 1993), with the remainder likely a product of some combination of SNe Type la 
and Type II. Both the weak and main s-process are "secondary" processes in the sense 
that a "secondary" element requires a Fe (or Fe-peak) seed nucleus to build upon. Such 
an element is expected to decrease in abundance (relative to Fe) as the star's overall iron 
content decreases (because the production of the element depends upon the availability of 
seed Fe nuclei). In contrast, a "primary" process does not require a Fe seed nucleus, and 
so the abundance of a primary element (relative to iron) should be constant with overall 
metallicity. While [Cu/Fe] does decrease with metallicity in the manner of a secondary 
element, the s-process contribution to copper is sharply restricted. Increasing the s-process 
component of copper in s-process nucleosynthesis models overproduces other light s-process 
elements (like Sr), and models of chemical evolution fail to reflect the observed solar 
composition (Matteucci et al. 1993; Baraffe & Takahashi 1993). The massive star (weak 
s-process) contribution appears to be an order of magnitude too small to reproduce the 
Cu content of either the Sun or metal-poor stars, and it must be augmented by explosive 
nucleosynthesis. It is unclear whether the additional copper in field stars is produced mainly 
by Type la supernovae (Matteucci et al. 1993; Baraffe & Takahashi 1993) or by Type II 
supernovae (Timmes et al. 1995). 

Cohen (1978, 1979, 1980) 's early observations of copper abundances in giant stars of 
several globular clusters of varying metallicity (—2.4 < [Fe/H] < —0.4) yielded somewhat 
ambiguous results. Cohen found evidence that copper is deficient with respect to iron 
and that this deficiency may become larger at lower metallicities. However, uncertainties 
in the derived abundances related to hyperfine broadening in the Cu line made a clearer 
delineation of any such trend impossible (Cohen 1980). Early studies such as these showed 
no real abundance trends in globular cluster iron-peak elements, copper included (e.g., 
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Spite & Spite 1985). 

The situation is further complicated by the behavior of [Cu/Fe] with respect to [Fe/H] 
in uj Centauri, a massive globular cluster that displays the largest confirmed range in stellar 
metallicity, with stars as metal poor as [Fe/H] ~ -1.8 and as metal rich as [Fe/H] ~ -0.8 
(Norris et al. 1996; Suntzeff & Kraft 1996). Not only is the [Cu/Fe ]-ratio of u Centauri 
stars uniformly low over its metallicity range, there is also no evidence of a decrease in 
[Cu/Fe] with decreasing [Fe/H] (Smith et al. 2000; Cunha et al. 2002) as is found for field 
stars over the same interval. Cunha et al. (2002) suggest that this could simply be a 
reflection of uj Cen's unique chemical history, and may indicate that the cluster built up 
copper through the contributions of Type II SNe only, contributions from Type la SNe 
being either suppressed or not retained by the cluster. 

In this paper we present Cu abundances from several other globular clusters spanning 
the metallicity range —2.4 < [Fe/H] < —0.8. Our purpose is to determine whether a large 
sample of globular cluster stars show the same trend in [Cu/Fe] that field stars do, or 
whether they follow a pattern similar to uj Cen. Cu abundances were derived for bright red 
giants from M3, M4, M5, M10, M13, M15, M71, NGC 288, NGC 362, and NGC 7006. A 
comparison of the [Cu/Fe] -ratios of the two clusters M4 and M5 is especially interesting. 
Although the overall metallicity of M4 is very similar to that of M5, it shows significantly 
higher abundances of the elements Si and Al, as well as overabundances of s-process 
elements like La and Ba compared to M5 (Ivans et al. 1999, 2001). Ivans et al. (2001) also 
determined that the abundances of M5 were "normal" for population II stars- that is, M5's 
[X/Fe] ratios are similar to those found in halo globular cluster and field stars. 
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2. Observations 

All spectra have been previously analyzed by the California/Texas group (hereafter 
"CTG"), and we adopt their published model parameters as summarized in Table 1. In 
these analyses, the raw spectra were trimmed, flat fielded, extracted, and wavelength 
calibrated according to standard procedures using the data analysis packages IRAF 1 and 
SPECTRE (Fitzpatrick & Sneden 1987). For details we refer the reader to the appropriate 
source (see Table 1). The data set is composed of contributions from several observatories 
and instruments, and the data quality and wavelength range are not uniform. Spectral 
resolution ranges from R~ 30,000 to R~ 60,000 and S/N is typically ~ 70 at the Cu I line 
at 5782 A. For lower S/N spectra, a two-pixel Gaussian smoothing was applied to the data. 

The two optical Cu I lines that are strong enough to be useful for abundance 
determinations in metal poor stars lie at 5782.14 A and 5105.50 A. No detectable ionized 
lines exist at optical wavelengths. The spectral region near the 5782 A transition is shown 
for representative stars from three clusters in Figure 1. This line is strongest in the most 
metal-rich cluster in the sample (M71) but essentially undetectable in the most metal-poor 
cluster (M15). Since wavelength coverage for our spectra is not uniform, the 5105 A line 
was available only for NGC 288, NGC 362, M4, and those M15 and M13 spectra taken 
with the Hamilton echelle spectrograph (at R~50,000). While it is preferable to have 
independent abundance determinations based on several different lines, the 5105 A line is 
likely to be a less reliable Cu abundance indicator in certain cases. It is an intrinsically 
stronger transition (log gf = —1.52; Bielski 1975) and so comes much closer to saturation 
as metallicity increases. This makes the derived abundance more uncertain, since significant 

TRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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changes in the Cu abundance produce only small changes in the appearance of the line. A 
strong line also gives an abundance that is more sensitive to the choice of microturbulence. 
In contrast, log gf = —1.72 for the 5782 A line, which arises from the same excitation level 
as the 5105 A line (Bielski 1975). 

Besides these intrinsic problems, the 5105 A line also lies in a more crowded part of 
the spectrum, especially in the cooler and more metal rich stars. Several MgH features 
exist near 5100 A and become quite strong in the cooler stars even at low metallicity, 
making continuum identification and reliable synthesis fits more difficult. It is only in the 
most metal poor cluster studied here, M15 ([Fe/H]~ —2.4), that the MgH lines are absent. 
Finally, in these bright but cool giants the spectral energy distribution favors a much higher 
flux at red wavelengths than at blue wavelengths. This results in poorer S/N near the 
5105 A line relative to the 5782 A line. 



3. Copper Line Analysis 

3.1. Line parameters 

Copper has significant hyperfine components and two stable isotopes: 63 Cu and 
65 Cu. 63 Cu is the more abundant isotope in the solar system, contributing 69% of the 
total solar system abundance. The hyperfine and isotopic structures of the copper lines 
have significant broadening effects, which become more important as the Cu line gets 
stronger. We therefore chose to derive abundances from spectrum synthesis, using oscillator 
strengths and hyperfine coefficients given in the Kurucz line list (Kurucz 1993). Deriving 
the Cu abundance from equivalent widths would require that this broadening be somehow 
accounted for (e.g., artificial addition of microturbulence) with attendant increases in 
abundance uncertainties (Cohen 1978). Additionally, in the crowded spectral regions of cool 
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or metal rich stars, a given line is likely to be blended to some extent, resulting in additional 
uncertainty in an EW analysis. We therefore used the current version of the LTE line 
analysis program MOOG (Sneden 1973) to fit a synthetic spectrum to a 10 A wavelength 
region around the 5782.14 A Cu I transition. For those clusters for which we had sufficient 
wavelength coverage, a 10 A spectrum synthesis was also employed to derive an abundance 
from the Cu I line at 5105.5 A. 

The model atmospheres and parameters used in the spectrum synthesis are identical to 
those derived by the California/Texas group. We have adopted those previously determined 
abundances, particularly Mg and C, since the 5105 A region contains a large number of 
MgH and C 2 lines in the cooler and higher metallicity stars. Though these molecular lines 
have been treated as a simple parameter of fit, in only one case has it been necessary 
to alter either of these abundances. In the M4 star L2208 the Mg abundance derived 
from atomic Mg lines by Ivans et al. (1999) did not seem to reproduce the strong MgH 
features surrounding the 5105 A line. In this case, an increase of 0.3 dex to Ivans et al.'s 
Mg abundance better fit the observed features. This star was problematic in the Ivans et 
al. study (see their §4.2.2 for a discussion of the anomalies found in this star's alpha and 
light odd-elements. The data for this star are of lower resolution and lower S/N than most 
of the sample). However, since the Mg abundance serves here only as a parameter of the 
spectrum synthesis fit, it does not affect the Cu abundance we derive. 

3.2. The Solar Cu Abundance 

A line list was prepared by using MOOG (Sneden 1973) to fit lines taken from the 
Kurucz (1993) line list in a 10A region around each Cu line to the Kurucz et al. (1984) 
solar spectrum. While the fit was in general quite good with the initial line parameters, gf 
values for a few Fe lines were modified empirically to produce a better fit. The Cu I gf 
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and hyperfine values are well determined and were not changed. Smoothing parameters and 
the continuum level were adjusted to minimize the difference between the observed solar 
spectrum and the computed spectrum. Once the best fit to the overall spectral region had 
been achieved, the abundance of Cu was varied to produce the best fit to the Cu I line, and 
a solar abundance was determined for both Cu I transitions. 

With the Holweger- Mueller solar model (v t = 1.15; Holweger & Mueller 1974) we 
derived loge G = 4.06 from the 5782 A line, which agrees well with the value Cunha et al. 
(2002) derive in a similar manner (loge & = 4.06), but is lower than the accepted solar 
photospheric abundance (loge Q = 4.21 ± 0.04; Anders & Grevesse 1989; Grevesse & Sauval 
1998). The 5105 A line gives a larger abundance in better agreement with the standard 
value: loge Q = 4.21. Because of the disagreement, [Cu/Fe] 5 i 5 is calculated using the 
5105 A solar value and [Cu/Fe] 578 2 with the 5782 A solar value. This is done with the 
expectation that this will eliminate any systematic errors between the atomic parameters of 
the lines, and [Cu/Fe] from both lines can be compared without further correction 2 . Indeed, 
[Cu/Fejsios and [Cu/Fe]s782 for cluster stars generally agree when calculated in this manner 
(see Figure 2). 



3.3. Error Analysis 

The Cu abundances we derive are sensitive to the adopted stellar parameters T e ^, 
log g, v t , and [Fe/H]. We find that a change in T e ff of ± 100 K produces the largest effect 
in [Cu/H], a change of =p 0.10 dex. [Cu/H] is less sensitive to the assumed gravity: a change 
in logg of ± 0.2 dex produces a change in [Cu/H] of ± 0.04 dex. However, these changes are 

2 As in other CTG work, we have adopted log e Q Fe = 7.52 throughout; see Sneden et al. 
1991b for a discussion of this choice. 
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largely canceled out in [Cu/Fe] by the accompanying change in [Fe/H] (see, for example, 
Ivans et al. 2001, Table 3). Factoring in the change in [Fe/H] with T e ff and log g results 
in A [Cu/Fe] ±0.02 dex for A T e jj± 100K and no discernible change in [Cu/Fe] with a 
0.2 dex change in log g. A change of ± 0.20 dex in [Fe/H] produces a change of ± 0.04 in 
[Cu/H]. The effect of microturbulence varies, being essentially negligible in warm stars to 
effecting changes of up to =f= 0.10 dex in [Cu/H] for Av t = ±0.20 km/s in the coolest stars. 
The quality of the synthesis fit to the observed spectrum also varies, depending on S/N. In 
high S/N spectra (e.g., M4 and M5), the fit is good to about 0.05 dex. In other spectra, 
the fit is good to less than 0.10 dex. We have therefore adopted an error of 0.10 dex on all 
our abundances, though in some cases the synthesis fit may be good to only about 0.20 dex 
(regardless of other error sources). Those stars with larger fitting uncertainties are marked 
in Table 1. Additionally, the weak lines in the cluster M15 (while advantageous in some 
respects) make abundance determinations more uncertain, and we give all those stars an 
error of ± 0.20 dex on the basis of spectrum fitting alone. 

4. CTG Analysis of Globular Cluster Giants 

4.1. Methodology 

Methods of deriving abundances in old metal-poor stars vary considerably from one 
group of investigators to another. In §6 we will compare values of [Cu/Fe] obtained here for 
globular cluster giants with those obtained by Mishenina et al. (2002) for halo field giants, 
so we state at this point the assumptions that have gone into the CTG modeling, and a 
corresponding discussion of the Mishenina et al. methods can be found in Appendix A. 

CTG proceeds as follows: 

(1) T e ff. A preliminary estimate is made from knowledge of the cluster color-magnitude 
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diagram and estimated color excess (e.g., Harris et al. 1999), coupled to a color vs. T e // 
scale based on (B-V) or preferably (V-K) (e.g., Cohen 1978, Alonso et al. 1999, Houdashelt 
et al. 2000). The preliminary value is modified with the requirement that the value of 
[Fe/H] derived from the LTE analysis of each Fe I line shall be independent of the line's 
lower excitation potential. The adopted value of T e ff is thus essentially independent of 
errors in color or color excess. 

(2) log g. It is assumed that LTE is satisfied and that logeiYe) from Fe II can be 
brought into agreement with logeiYe) from Fe I. This technique avoids uncertainties in the 
distance scale associated with globular clusters (see, e.g., Reid 1997, 1998 and Carretta 
et al. 2000). However, Thevenin & Idiart (1999) have pointed out that NLTE effects in 
stellar atmospheres lead to the over-ionization of iron, such that measurements of Fe I 
and Fe II need not result in the same Fe abundance. This effect has implications for the 
derivation of both log g and [Fe/H]. 

Since both [Cu/H] and [Cu/Fe] are insensitive to uncertainties in log g (see §3.3), 
an LTE derivation of this quantity will not affect those results. The remaining issue, the 
metallicity scale for these clusters, is beyond the scope of this paper but has been addressed 
by Kraft & Ivans (2002). Kraft & Ivans (2002) follow the reasoning of Thevenin & Idiart 
(1999) and Ivans et al. (2001), who argue that Fe II is a better measure of stellar iron 
abundances than Fe I (because most of the iron atoms in the star are ionized), and measure 
globular cluster metallicities accordingly. We have preserved the original LTE metallicities 
as measured by CTG, under the assumption that the effects of over-ionization on Cu I in 
a star's atmosphere are somewhat mitigated if one measures abundance ratios of elements 
with similar ionization potentials in the same ionization state (i.e., deriving [Cu/Fe]/; see 
Kraft & Ivans 2002 for a discussion of these issues). In an LTE analysis, the reported 
metallicity is essentially the Fe I abundance and our [Fe/H] values differ from Kraft & 
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Ivans (2002) by varying amounts. However, maintaining the LTE [Fe/H] scale preserves our 
ability to directly compare the run of [Cu/Fe]/ with [Fe/H] in clusters with that of the field 
and uj Cen, since those studies have also employed the standard LTE analysis. NLTE effects 
will distort the [Fe/H] axis of our trends, but legitimate comparisons to field star trends 
can still be made so long as we take care to match cluster giant stars with field giant stars. 

(3) Values of log gf. For Fe I, the adopted values are consistent with the Oxford group 
(Blackwell et al. 1982) and with Gurtovenko & Kostyk (1989) (see also Fuhr et al. 1988). 
For Fe II, CTG adopted those of Moity (1983), Heise & Kock (1990) and Biemont et al. 
(1991). A more nearly complete discussion is given in Sneden et al. (1991a). In M4 and 
M5, these lines were supplemented by revised log g/'s from Lambert et al. (1996). 

(4) Models. CTG employs MARCS models for old, metal-poor giants (Gustafsson et 
al., 1975). In connection with this exposition of method, we make a special comment on 
the CTG analysis of giants in M5 and M4. In the M5 paper Ivans et al. (2001) prefer 
to adopt a procedure in which T e ff is derived from reddening corrected colors and the 
Alonso et al. (V — K) vs. Tg/j scale, and log g is derived from knowledge of the cluster 
distance modulus, this in response to concerns about the reliability of the assumption of 
LTE (compare Thevenin & Idiart 1999). This alters somewhat the derived abundance 
of Fe, in comparison with the (above) "traditional" method. This approach is, however, 
inconsistent with the CTG analysis of the other clusters. Fortunately, Ivans et al. also give 
the result of the traditional analysis (see Appendix A, Table 11), and we adopt here the 
value of [Fe/H] ai , 9 = —1.35 given in the appendix. For M4, which is discussed in the M5 
paper in connection with this problem, we also adopt here the metallicity derived from the 
"traditional" approach: [Fe/H] ai , 9 = —1.18. 

Model parameters, values of [Fe/H] taken from CTG, and values of [Cu/Fe] derived 
here for individual cluster stars in M71, M4, M5, NGC 362, NGC 288, M3, M10, NGC 7006, 



- 12 - 



M13 and M15 are given in Table I. Cu abundances are normalized to the mean of [Fe/H] 
derived from Fe I and Fe II. Some clusters give rise to specific concerns, and we address 
those in the following section. 

4.2. Special considerations 

4.2.1. M4 

M4 was one of the few clusters in our sample for which both Cu lines were available for 
analysis. However, the 5782 A feature falls close to the end of the spectral orders in this 
cluster, and so ~ 4A of the spectrum used for fitting in the other clusters is unavailable. 
This is not likely to be a source of significant error since the spectrum in these stars is very 
clean (high S/N and relatively weak features). 

A potentially more serious problem in M4 is the presence of a diffuse interstellar band 
(DIB) at 5780A (e.g. Herbig 1975; Krelowski & Sneden 1993). Fortunately, the Cu feature 
is shifted enough by M4's radial velocity that the DIB does not directly affect it (see Figure 
3). To improve the fit of the synthetic spectrum, the DIB was divided out manually using 
SPECTRE (Fitzpatrick & Sneden 1987), effectively treating the feature as continuum. 
Comparison of [Cu/Fe] 578 2 to [Cu/Fe] 5 i 5 indicates that the DIB did not in fact affect the 
Cu feature and its removal by this simple means has not introduced any appreciable fitting 
error; the spectrum fits for this line are as good as ±0.05 dex. In every case but two (stars 
L2208 and L4201), the abundances from the 5782 A and 5105 A lines agree to 0.1 dex. 

Though the 5105 A line was available for all M4 stars, it was not possible in every 
case to determine an abundance from it. In stars cooler than 4225 K (with the exception of 
L2406) the spectrum was too crowded by strong MgH lines to achieve a reliable fit of the 
synthesis to the observed spectrum. 



4-2.2. M15 



M15 is the most metal poor cluster in this sample. The overall metallicity is so low 
that the 5782 A feature is too weak to synthesize even in the cooler stars; all that could be 
derived from it were unenlightening upper limits. We have therefore synthesized only the 
5105 A line in M15 stars. Given the usually good agreement between the two lines in other 
clusters, we are confident that the derived abundance from this line accurately reflects the 
Cu content of a particular star. However, even the stronger 5105 A line is still quite weak 
in the M15 spectra, and individual abundances are still uncertain. We put larger errors on 
M15 abundances, and our average abundance for this cluster has a high a. 



4.2.3. Mil 

The most metal rich cluster in the sample, M71, posed a different sort of problem 
for our analysis. The spectrum becomes so crowded with lines that it is very difficult to 
determine the level of the continuum. Only six stars from Sneden et al. (1994) yielded 
credible results. In the coolest stars, we found that [Cu/Fe] correlated with T e //, with the 
coolest stars giving the lowest [Cu/Fe]. This is likely due to errors in the continuum fitting. 
In the coolest giants, molecular band heads of TiO are particularly strong (see Sneden et al. 
1994, their Fig. 1), and a TiO band head does exist near the 5782 A feature (at 5758 A, 
degraded to the red; Pearse & Gaydon 1963). If these band heads are not properly taken 
into account, the tendency is to incorrectly identify the continuum so that the Cu line 
becomes shallower than it ought to be and gives a falsely low Cu abundance. 
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5. [Cu/Fe] Ratios in Globular Clusters 

Copper lies at the end of the Fe-peak and may or may not actually be formed in the 
same processes that create the lighter elements of the Fe-peak. However, like the Fe-peak 
elements, it shows no real star to star variation in these globular clusters. The standard 
deviation of the mean [Cu/Fe] value for every cluster in the sample is quite small and well 
within the observational errors. In Table 2 we list for each cluster the values [Fe/H] a „ s , 
[Cu/Fe] a „ s , and the standard deviation (a), and in Figure 4, we plot [Cu/Fe] a „ 9 vs. [Fe/H] a „ s 
for the clusters in our sample. 

As a supplement to Figure 4, we offer in Figure 5 box plots in which the clusters are 
arranged in order of increasing metallicity, and which illustrate both the median [Cu/Fe] 
and the range in [Cu/Fe] for each cluster. Figure 5 appears to show a rather large range in 
[Cu/Fe] in M10, when the outliers are considered. However, an examination of the spread 
in [Ni/Fe] (which ought not to vary from star to star) in M10 as measured by Kraft et al. 
(1995) indicates that the spread in [Cu/Fe] is not real: an investigation of each of the five 
Ni lines analyzed by Kraft et al. (1995) reveals that cr( [Ni/Fe]) ranges from 0.10 dex to 
0.15 dex. <r([Cu/Fe])=0.10 is therefore not significant. The total range of [Ni/Fe] in each 
individual Ni line is as much as 0.3 dex in this cluster. It should also be noted that the 
single low outlying [Cu/Fe] value, M10:G, is the warmest star in our M10 sample, and has 
a very weak line and low S/N. Kraft et al. (1995) were able to measure only one Ni line 
in that star. In M4, a cluster for which we have a large number of high quality spectra, 
cr( [Ni/Fe]) is also comparable to <r([Cu/Fe]), which suggests that the spread in [Cu/Fe] in a 
particular cluster is related to data quality rather than any actual star to star variation. 

Among halo field stars, which are generally regarded as surrogates of globular cluster 
stars, Sneden et al. (1991a) found a sharp decrease in [Cu/Fe] with decreasing metallicity 
(a mean of Fe-peak elements including Ni but dominated by Fe), following the mean 
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relationship [Cu/M]= +0.38[M/if] + 0.15. We plot this also in Figure 4. It is seen that in 
general the clusters follow a fairly similar decline with roughly the same slope as that of 
the field stars, but displaced toward slightly higher values of [Fe/H] and/or lower values of 
[Cu/Fe]. With the exception of M15, all clusters lie "below" the line given by Sneden et al. 
(1991a). 

The displacement between field and cluster stars noted from inspection of Figure 4 may 
possibly be real, but may also be a result of systematic differences between investigators. 
We explore this possibility further using the expanded sample of halo field stars recently 
studied by Mishenina et al. (2002), who explored [Cu/Fe]-ratios in a sample of 90 metal-poor 
stars in the range —3.0 < [Fe/H] < —0.5, on the basis of high-resolution, high S/N spectra. 

6. Clusters vs. Field 

Out of the 90-star sample of Mishenina et al., we first selected for comparison with 
cluster giants the 21 low metallicity ([Fe/H] < —0.50) field giants having log g < 2.3. 
These are essentially the field star surrogates of the cluster giants; elimination of the others 
suppresses any possible systematic effects owing to a large shift in log g and those NLTE 
effects discussed in §4.1. In addition, the lower luminosity sample is heavily weighted 
toward stars having metallicities [Fe/H] > —1.0, a regime more nearly associated with the 
thick disk than with our (predominantly) halo globular clusters. 

Before making the comparison of [Cu/Fe] values among these 21 giants with our 
globular cluster giants, we must first eliminate as far as possible systematic effects. The 
details of their abundance method and the resulting offsets from this work are described in 
Appendix A, and those results are listed in Table 2. 

In Figure 6 we plot [Cu/Fe] vs. [Fe/H] for both clusters and halo field giants, having 
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adjusted the clusters to system of abundances defined by Mishenina et al. (Obviously, we 
could have normalized equally well to the system defined by CTG.) The shifts required to 
move all stars to the same abundance system are small: this is largely because the methods 
of analysis are quite similar. Inspection of Fig 6 shows that the cluster and halo field giants 
have a nearly identical distribution of [Cu/Fe] as a function of [Fe/H]. From [Fe/H] near 
—0.7, the value of [Cu/Fe] is near —0.1, then it declines to about —0.65 at [Fe/H]~ —1.7, 
after which it levels off at a constant value near —0.65 as [Fe/H] approaches —3.0. 



7. Conclusion 
7.1. M4 and M5 

The two clusters M4 and M5, while very similar in [Fe/H], show a significant difference 
in their copper abundance. This is not likely to be a result of error in the analysis or in 
stellar atmospheric parameters, given the very small range in [Cu/Fe] derived for these 
clusters over a large range in T e ff and log g. Furthermore, the difference in the copper 
abundance between similar stars from different clusters reinforces the difference seen on 
average. Figure 7 illustrates this with L2206 and IV-34, stars typical of the M4 and M5 data 
sets. These two stars have very similar atmospheric parameters, yet a simple examination of 
their spectra shows that Cu is indeed overabundant in M4. While (with one exception) no 
other species in this region shows significant line depth contrasts, the flux at the 5782 A line 
in the M4 star is ~1.3 times weaker than in the M5 star. Comparing this line with others, 
it is evident that the effect is much greater than simple noise or continuum uncertainties. 
The other excess present in Figure 7 is Si, an element shown by Ivans et al. (1999) to be 
overabundant in M4. Unlike other elements in M4, however, Cu does not appear to be 
unusually enhanced compared to field halo stars. In this metallicity region, the field halo 
stars show a very swift increase in [Cu/Fe] that appears to be mirrored in M4, and continues 
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in the higher metallicity clusters M71 just as it does in the field. Whatever M4's source of 
enrichment may be that sets it apart in other elements, it has not significantly contributed 
to the Cu abundance (any more than is "normal" for its [Fe/H]). 

Like M4/M5, NGC 288/NGC 362 have similar metallicities. NGC 288 and M4 have 
similar Na-0 and Al-0 anticorrelations while M5 and NGC 362 have similar Na-0 and Al-0 
anticorrelations (Shetrone & Keane 2000). Unlike the M4/M5 pair, the NGC 288/NGC 32 
pair does not show such a large discrepancy in [Cu/Fe]. Though NGC 288 is elevated in 
Cu, this effect is smaller than in M4, and the uncertainty in NGC 288's Cu abundance 
(er = 0.11) is large. Though NGC 362's Cu abundance is relatively well determined, that of 
NGC 288 is not, and we can draw no conclusions as to Cu abundance differences between 
the two. These two clusters also seem to follow the field halo star Cu abundance trend, 
though M4, M5, NGC 288, and NGC 362 all fall near the [Fe/H] value where the trend in 
[Cu/Fe] with [Fe/H] changes shape in the field (see Fig. 6). 

7.2. Clusters vs. lu Cen 

Cluster and field halo giants both display a trend in [Cu/Fe] with [Fe/H] that is totally 
absent in u Cen, where [Cu/Fe] is uniformly deficient over a wide range in [Fe/H]. At or 
about the metallicity of M5 (~ —1.35), [Cu/Fe] begins to increase very quickly with [Fe/H] 
in the field and in the higher metallicity clusters. 

Ivans et al. (1999) found overabundances of s-process elements like Ba and La in 
M4. uj Cen also shows evidence of s-process enhancement in its more metal-rich stars 
(Cunha et al. 2002). It would be tempting to attribute M4's increased Cu abundance 
(and therefore also that of the field) to enhanced s-processing as well. However, the bulk 
of solar system Cu is made from the weak, not main, s-process, unlike Ba and La. Only 
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about 10% of the solar system Cu can be attributed to the main s-process (Kappeler et al. 
1989). Baraffe & Takahashi (1993) find that the massive star (weak s-process) contribution 
is an order of magnitude too small (for all [Fe/H]) to reproduce the field star abundance, 
and increasing the main s-process contribution to M4 should not make much difference in 
the Cu abundance. In any case, even s-process enhanced stars in cu Cen still show a low 
[Cu/Fe] (though there is some Cu enhancement of giants at the highest metallicities; see 
Pancino et al. 2002). Cunha et al. (2002) attribute the low [Cu/Fe] in cu Cen to enhanced 
SNe Type II (and little SNe Type la) enrichment. If so, then the enhancement in M4's Cu 
cannot be due to Type II SNe, either. This conclusion is reinforced by the fact that the 
Cu abundances in M15 stars do not correlate with the Eu abundances. Eu, an r-process 
element, is commonly associated with Type II SNe. If Cu owed its production to the 
same process, then it ought to show the same star-to-star variation seen in the M15 Eu 
abundances. An inspection of the Eu abundances in Sneden et al. (1997) for M15 stars 
shows that the highest and lowest [Eu/Fe] stars in our sample, K462 and K583, have 
virtually identical [Cu/Fe]. The highest and lowest [Cu/Fe] stars, K634 and K853, have 
essentially the same [Eu/Fe]. There is therefore no indication that [Cu/Fe] is in any way 
related to the SNe Type II tracer [Eu/Fe] in M15. 

Taken together, the Cu abundances in globular clusters and halo field stars suggest 
a third possibility, that at higher metallicity stellar Cu abundances receive a significant 
contribution from Type la SNe. The abundance pattern of cu Cen, where [Cu/Fe] is low, 
is likely the product of Type II SNe (Smith et al. 2000; Cunha et al. 2002), making it 
improbable that Cu comes from this source, as Timmes et al. (1995) suggest (though 
Timmes et al. do show that the copper yield of SNell is metallicity-dependent). Substantial 
s-process contributions are ruled out on the basis that M4, a cluster known to be s-process 
enhanced, shows no evidence for an unusual Cu enhancement. Matteucci et al. (1993) have 
suggested that Cu may be produced by SNe la, a source that seems increasingly likely. 
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Possibly the best way to resolve this question is to further study M4 and M5's abundance 
patterns; these two clusters lie at a metallicity where the nucleosynthetic source of Cu 
appears to undergo a transition. In particular, the relative abundances of Zn (the next 
heavier element) in these two clusters may shed some light on the origin of Cu. 
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A. Appendix 

In order to determine systematic offsets between our cluster analysis and the field star 
analysis by Mishenina et al. (2002), we take up the procedures employed by Mishenina et 
al. and invite the reader to compare these with the CTG methods described in §4.1. The 
actual stellar models, however, are not discussed in Mishenina et al. (2002), but rather are 
adopted from an earlier paper, Mishenina & Kovtyukh (2001). 
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(1) Equivalent widths. These are in excellent agreement with those obtained by 
Pilachowski et al. (1996) (see Mishenina & Kovtyukh 2001, their Fig 3), which are on the 
system of CTG. 

(2) T e ff- The models are assigned preliminary values of T e ff from colors and the 
Alonso et al. scale of color vs. T e ff- However, final assignment of T e ff depends upon 
achieving "flatness" in the Fe I excitation plot. The procedure for setting T e // in the end is 
identical to that of CTG. 

(3) log g. The abundance of Fe from Fe II is brought into agreement with that of Fe 
from Fe I to set log g, again a procedure identical with that adopted by CTG. 

(4) log gf values. For Fe I and Fe II, Mishenina et al. adopt the log gf values (as 
do Mishenina & Kovtyukh) of Gurtovenko & Kostyk (1989). For the 30 Fe I lines in 
common between CTG and Mishenina & Kovtyukh, we find the remarkable result that 
A log gf = 0.000 (a = 0.066). For the six Fe II lines employed by CTG, there is a small 
but significant offset, A log gf = +0.108 (a = 0.064), in the sense CTG minus Mishenina 
& Kovtyukh. We may therefore expect to find a difference in the values of log g between 
the present investigation and that of Mishenina et al. 

For Cu, Mishenina et al. adopted the log gf values given by Gurtovenko & Kostyk 
(1989) from inverted solar analysis. For the Cu line at 5782 A, their log gf value is identical 
with that of CTG. For the line at 5105 A, their adopted log gf value is smaller than that 
of CTG by 0.06 dex. Thus from this line alone, they would derive a larger abundance of Cu 
than we derive in the present investigation; this amounts to 0.03 dex, averaging over the 
two lines. Although we do not apply this "correction" in the comparison of results, we note 
that 5105 A is analyzed in this paper only for a few giants in M4, NGC 362, NGC 288, M13 
and all stars in M15. 
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(5) Adopted models. Mishenina et al., following Mishenina & Kovtyukh, adopt 
Kurucz models with convective overshoot turned on (Kurucz 1993). The differences 
induced in [Fe/H] among metal-poor giants by switching from MARCS to Kurucz models 
has been explored by Kraft & Ivans (2002). Among intermediately metal-poor stars 
(-1.7 < [Fe/H] < - 1.2), the effect on [Fe/H] (from Fe I) is negligible (< 0.01 dex), but 
rises to 0.07 to 0.10 dex among relatively metal-rich and metal-poor giants. 

Fortunately, of the 21 halo field giants considered here, 11 have also been studied by 
CTG (Kraft et al. 1992; Shetrone 1996). We can then empirically determine the offsets in 
T e jj, log g, and [Fe/H] which must be applied to bring the field and cluster giants onto 
the same abundance system. In the sense CTG minus Mishenina & Kovtyukh, we find the 
following mean differences: A T e// = +A7K (a = 89.5/T); A log g = +0.28 (a = 0.20); 
A [Fe/H]= +0.05 (a = 0.10). 

The offsets in T e // and [Fe/H] are small (if significant), but the offset in log g is rather 
large. Fortunately, the values of [Cu/Fe] are not very sensitive to variations in log g. In 
addition, and more importantly, the offset is just about what one would expect owing to 
the differences in the log gf values for Fe II. If [Fe/H] from Fe II is brought into agreement 
with [Fe/H] from Fe I, then a shift of 0.11 in the Fe II log gf values corresponds to a shift 
in log g of 0.20 dex (Ivans et al. 2001, Table 3), close to the value of 0.28 cited above. The 
[Fe/H] offset cited above is also explicable in terms of the offset in T e ff- Again, according 
to the Ivans et al. Table 3, an increase in T e ff of 47K should generate an increase of 0.04 
dex in [Fe/H] from Fe I, almost exactly that observed. 

Having made the appropriate adjustments in [Fe/H] and T e ff, it remains only to adjust 
the values of [Fe/H] from one set of atmospheric models to the other. In columns 5 and 6 
of Table 2, we list also mean cluster values of [Fe/H] and [Cu/Fe] on the basis of Kurucz 
models. There is almost no effect for most of the clusters, but there are significant shifts in 



the case of M71 and M15. 
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Fig. 1. — Representative Cu I 5782 A lines in stars of three clusters: M71 (upper panel), 
M13 (middle panel), and M15 (bottom panel). The 5782.14 A line is marked in each. M71 
and M15 are the metallicity extremes in the sample. Most of the clusters have [Fe/H] similar 
to M13. 
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Fig. 2. — Comparison of copper abundances derived from the 5782A and 5105A lines. In 
most cases abundances from the two lines agree to within 0.10 dex. 
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Fig. 3.— Smoothed spectrum of M4:L2206 (dotted line) and M5:IV-34 (solid line). The stars 
have similar T e ff, log g, [Fe/H], and microturbulence. The dip in the M4 star's continuum 
is due to a DIB near 5780 A, though it clearly does not extend into the Cu I feature. 
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Fig. 4. — Average cluster [Cu/Fe] values on the MARCS model scale. The line indicates 
the trend found by Sneden et al. (1991a). Most clusters fall below the trend, but the overall 
shape appears similar. 
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Fig. 5. — The boxed region for each cluster encompasses the interquartile range (middle 
50%) of the [Cu/Fe] data for that cluster. Also shown are the median (horizontal line), 
range (vertical lines; excludes outliers) and outliers (crosses; an outlier has a value greater 
than 1.5 time the interquartile range). The clusters are arranged in order of increasing 
[Fe/H] ai , 9 . These values are not corrected to the Mishenina et al. system. 
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Fig. 6. — Cluster averaged [Cu/Fe] corrected onto the system used by Mishenina et al. 
(2002). The plotted field stars are giants, in keeping with the CTG sample. There is no 
indication that the cluster stars deviate from the field field stars in [Cu/Fe] at a comparable 
[Fe/H]. This also holds true for M4, a cluster that is overabundant in Ba, La, and Si. 
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Fig. 7. — Upper panel: Smoothed, interpolated spectrum of M4:L2206 (dotted line) and 
M5:IV-34 (solid line), two stars with similar T e ff, log g, [Fe/H], and microturbulence. Lower 
panel: the division of the M5 star by the M4 star, illustrating the excess depth of the M4 
star's Cu I line. The positions of two Si I lines have also been marked, affirming Ivans et 
al.'s claim of an excess of Si in M4. 
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0.75 


1.60 


-1.33 




-0.58 


-0.58 


IV-34 


4275 


0.65 


1.55 


-1.32 




-0.58 


-0.58 


1-58 


4350 


0.80 


1.50 


-1.31 




-0.56 


-0.56 


1-61 


4400 


1.00 


1.50 


-1.30 




-0.57 


-0.57 


11-59 


4463 


1.15 


1.65 


-1.25 




-0.61 a 


-0.61 a 


111-18 


4475 


0.55 


1.70 


-1.46 




-0.59 


-0.59 


1-2 


4500 


1.10 


1.45 


-1.36 




-0.57 


-0.57 


1-50 


4525 


1.15 


1.40 


-1.35 




-0.54 


-0.54 


11-74 


4525 


1.30 


1.30 


-1.25 




-0.67 


-0.67 


111-59 


4575 


1.20 


1.35 


-1.31 




-0.53 


-0.53 


IV-36 


4575 


1.50 


1.35 


-1.30 




-0.63 


-0.63 


11-50 


4590 


1.57 


1.45 


-1.20 




-0.66 


-0.66 


IV-30 


4625 


1.00 


1.75 


-1.44 




-0.63 


-0.63 


111-52 


4625 


1.50 


1.45 


-1.34 




-0.54 


-0.54 


IV-4 


4625 


1.55 


1.20 


-1.25 




-0.70 a 


-0.70 a 


IV-26 


4650 


1.05 


1.40 


-1.40 




-0.55 


-0.55 


111-53 


4700 


1.05 


1.75 


-1.46 




-0.63 


-0.63 
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Star T eff logg v t (km/s) [Fc/H] [Cu/Fe] - [Cu/Fe] - [Cu/Fe] at , 9 



1-55 


4700 


0.85 


1.80 


-1.48 




-0.57 


-0.57 


NIGC 288 (Shctrone & 


Keane 2000) 










531 


3780 


0.10 


1.60 


-1.31 




-0.50 a 


-0.50 


403 


3950 


0.20 


1.90 


-1.43 




-0.37 


-0.37 


274 


4025 


0.70 


1.90 


-1.37 




-0.33 


-0.33 


344 


4180 


0.80 


1.60 


-1.36 




-0.68 


-0.68 


245 


4250 


0.80 


1.40 


-1.41 


-0.10 


-0.40 


-0.40 


231 


4300 


1.10 


1.50 


-1.41 


-0.25 


-0.32 


-0.29 


388 


4325 


1.30 


1.60 


-1.37 




-0.32 


-0.32 


297 


4330 


1.20 


1.70 


-1.41 


-0.30 


-0.32 


-0.31 


351 


4330 


1.20 


1.55 


-1.33 


-0.30 


-0.48 


-0.39 


287 


4350 


1.20 


1.40 


-1.45 




-0.32 


-0.32 


307 


4350 


1.20 


1.35 


-1.40 




-0.45 


-0.45 


\13 (Kraft ct 


al. 1992) 














VZ-1397 


3950 


0.40 


2.35 


-1.46 




-0.50 


-0.50 


AA 


4000 


0.40 


2.25 


-1.45 




-0.55 


-0.55 


11-46 


4000 


0.60 


2.10 


-1.46 




-0.65 


-0.65 


VZ-297 


4070 


0.70 


2.25 


-1.52 




-0.55 


-0.55 


111-28 


4160 


0.75 


1.75 


-1.49 




-0.65 


-0.65 


VZ-1000 


4175 


0.45 


2.10 


-1.45 




-0.50 


-0.50 


VZ-1127 


4225 


0.90 


2.00 


-1.48 




-0.65 


-0.65 


V110 (Kraft ct al. 1995) 












A-II-24 


4050 


0.10 


2.00 


-1.50 




-0.69 


-0.69 


A-III-21 


4060 


0.50 


2.10 


-1.49 




-0.50 


-0.50 


H-I-367 


4135 


0.60 


1.70 


-1.50 




-0.72 


-0.72 


A-III-16 


4150 


0.90 


2.00 


-1.50 




-0.47 a 


-0.47 


B 


4150 


0.50 


1.80 


-1.49 




-0.64 


-0.64 


C 


4200 


0.75 


2.00 


-1.62 




-0.60 


-0.60 


D 


4200 


1.05 


2.00 


-1.49 




-0.55 


-0.55 


H-I-15 


4225 


0.75 


1.75 


-1.50 




-0.65 


-0.65 


E 


4350 


0.80 


2.00 


-1.60 




-0.70 


-0.70 


A-I-60 


4400 


1.10 


1.60 


-1.50 




-0.65 


-0.65 
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Star 




log g 


vt(km/s) 


[Fe/H] 


[ Cu /Fe] 5105A 


t Cu /Fc] 5782A 


[Cu/Fe] a „g 


A-III-5 


4400 


1.20 


1. to 


1 Qfi 




n fin 
-U.DU 


n fin 

-U.DU 


A-I-61 


4550 


1.00 








n 7n a 


n 7n a 

-u. 1 u 


G 


4650 


1.20 


1 R^i 
i-.OO 


1 ^fi 




-u.oo 


-u.oo 


NGC 7006 


(Kraft et al. 1998) 












1-1 


3900 


0.10 


2.25 


-1.50 




-0.65 


-0.65 


V19 


4100 


0.3 


2.4 


-1.57 




-0.65 


-0.65 


11-103 


4200 


0.75 


1.85 


-1.54 




-0.50 


-0.50 


11-46 


4200 


0.5 


1.75 


-1.56 




-0.55 


-0.55 


M13 (Kraft ct al. 


1997) 












L-598 


3900 


0.00 


2.10 


-1.62 


-0.57 


-0.55 


-0.56 


1-48 


3920 


0.30 


2.00 


-1.57 


-0.60 


-0.67 


-0.64 


L-629 


3950 


0.20 


2.00 


-1.64 


-0.75 


-0.58 


-0.67 


11-67 


3950 


0.20 


2.10 


-1.52 


-0.67 


-0.67 


-0.67 


11-90 


4000 


0.30 


2.00 


-1.58 


-0.50 


-0.65 


-0.58 


IV-25 


4000 


0.15 


2.25 


-1.56 


-0.45 


-0.75 


-0.60 


L-835 


4090 


0.55 


1.90 


-1.52 




-0.65 


-0.65 


1-12 


4600 


1.50 


1.60 


-1.58 




-0.60 


-0.60 


IV- 19 


4650 


1.50 


1.60 


-1.60 




-0.60 


-0.60 


11-09 


4700 


1.70 


1.50 


-1.58 




-0.58 


-0.58 


IV-22 


4700 


1.90 


1.50 


-1.54 




-0.60 


-0.60 


11-41 


4750 


2.00 


1.75 


-1.48 




-0.65 


-0.65 


1-72 


4850 


1.90 


1.40 


-1.64 




-0.40 b 


-0.40 b 


11-01 


4850 


2.10 


1.25 


-1.56 




-0.60 


-0.60 


11-28 


4850 


1.75 


2.00 


-1.64 




-0.65 a 


-0.65 a 


1-54 


4975 


1.70 


1.75 


-1.68 




-0.55 b 


-0.55 b 


M15 (Sncden ct al. 1997) 












K386 


4200 


0.15 


1.85 


-2.38 


-0.65 




-0.65 


K462 


4225 


0.30 


1.85 


-2.43 


-0.70 




-0.70 


K634 


4225 


0.30 


1.85 


-2.42 


-0.57 




-0.57 


K583 


4275 


0.30 


1.90 


-2.42 


-0.66 




-0.66 


K853 


4275 


0.50 


1.85 


-2.38 


-0.80 




-0.80 


K702 


4325 


0.25 


1.90 


-2.42 


-0.73 




-0.73 
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Star T eff logg v t (km/s) [Fc/H] [Cu/Fe] 5iosA [Cu/Fc] 57g2A [Cu/Fe] „ 

K479 4325 0.45 2.30 -2.42 -0.70 b ■•• -0.70 b 

K490 4350 0.60 1.65 -2.44 -0.45 ■•• -0.45 



a error =±0.2 dex 
b upper limit only 
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Table 2. Cluster average [Cu/Fe] 



Cluster [Fe/H] ovs [Cu/Fc] a „ 9 a [Ee/H] ot , g [Cu/Fc] a „ 9 [Fe/H] OU9 [Cu/Fe] avg 
MARCS MARCS Kurucz Kurucz corrected corrected 



M71 


-0.79 


-0 


17 


0. 


05 


-0.72 


-0 


.24 


-0.77 


-0. 


19 


M4 


-1.18 


-0 


32 


0. 


06 


-1.18 


-0 


32 


-1.23 


-0. 


27 


NGC 362 


-1.33 


-0 


52 


0. 


06 


-1.33 


-0. 


.52 


-1.38 


-0. 


.47 


M5 


-1.35 


-0 


.58 


0. 


08 


-1.35 


-0. 


58 


-1.40 


-0. 


53 


NGC 288 


-1.39 


-0. 


40 


0. 


11 


-1.39 


-0 


.40 


-1.44 


-0. 


35 


M3 


-1.46 


-0 


59 


0. 


07 


-1.45 


-0 


60 


-1.50 


-0. 


55 


M10 


-1.52 


-0 


64 


0. 


10 


-1.51 


-0 


.65 


-1.56 


-0. 


60 


NGC 7006 


-1.55 


-0. 


.59 


0. 


08 


-1.54 


-0. 


60 


-1.59 


-0. 


55 


M13 


-1.58 


-0 


.60 


0. 


06 


-1.57 


-0 


.61 


-1.62 


-0. 


56 


M15 


-2.41 


-0. 


66 


0. 


11 


-2.31 


-0 


76 


-2.36 


-0. 


71 



